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Abstract. We study the Deffuant et al model for continuous-opinion dynamics 
under the influence of noise. In the original version of this model, individuals meet 
in random pairwise encounters after which they compromise or not depending of a 
confidence parameter. Free will is introduced in the form of noisy perturbations: 
individuals are given the opportunity to change their opinion, with a given probability, 
to a randomly selected opinion inside the whole opinion space. We derive the master 
equation of this process. One of the main effects of noise is to induce an order-disorder 
transition. In the disordered state the opinion distribution tends to be uniform, while 
for the ordered state a set of well defined opinion groups are formed, although with some 
opinion spread inside them. Using a linear stability analysis we can derive approximate 
conditions for the transition between opinion groups and the disordered state. The 
master equation analysis is compared with direct Monte- Carlo simulations. We find 
that the master equation and the Monte-Carlo simulations do not always agree due to 
finite-size induced fiuctuations that we analyze in some detail. 

Keywords: Collective phenomena in economic and social systems: Interacting agent 
models. Non-equilibrium processes: Stochastic particle dynamics (Theory). 
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1. Introduction 

The application of techniques and tools from nonlinear and statistical physics to 
understand the dynamics of opinion changes in a society has become a topic of interest 
in recent years [Ij. A society can be thought of as a complex system composed by a 
large number of interacting individuals with diverse opinions. These opinions are not 
necessarily static, but they evolve due to a variety of internal as well as external factors, 
such as the influence of advertising and acquaintances, amongst others. As a result of 
this evolution, a consensus opinion could emerge (a vast majority on individuals adopting 
a similar opinion), or the population could fragment into a number of groups. To analyze 
the process of opinion formation, several models inspired from statistical mechanics have 
been developed. In those, the opinion held by an individual is a dynamical variable which 
evolves by some rules, usually with an important stochastic ingredient [2]. Models can 
be divided in two broad categories: discrete models where the opinion can only adopt a 
finite set of values [HJ HI [5j , and continuous models where the opinion of an individual is 
expressed as a real number in a finite interval [SI |71 El [9j . Discrete models are useful when 
analyzing cases in which individuals are confronted with a limited number of options 
(a political election, for example) where one is forced to choose amongst a finite set of 
parties. Continuous models are more suitable to analyze cases in which a single issue 
(legalizing abortion, for example) is being considered and opinions can vary continuously 
from "completely against" to "in complete agreement" . 

A continuous model introduced by Deffuant and collaborators [6] has received much 
attention recently. This model implements the bounded confidence mechanism by which 
two individuals only infiuence the opinion of each other if their respective opinions differ 
less than some given amount. In other words, people holding too distant opinions on 
an issue will simply ignore each other and will, hence, keep their original opinions. It 
is only through the interaction of not too distant people that we manage to modify our 
opinion. This model was in turn inspired by the Axelrod model for the dissemination of 
culture |10| and general threshold models [H] and has in turn inspired a large number 
of extensions and modifications [l2l[l3lEl[l5l[l6l[I7l[l8l[l9l[^ 

In the Deffuant et al. model individuals meet in random pairwise encounters in a 
given connectivity network, but the subsequent evolution is completely deterministic. 
This leads to final states in which either perfect consensus has been reached or the 
population splits in a finite number of groups such that all individuals in one group have 
exactly the same opinion. We believe that such uniform states are not very realistic and 
some degree of discrepancy must appear within otherwise well defined groups. In this 
paper, we introduce an additional element of randomness in the dynamics. It aims to 
represent, certainly in a caricaturist manner, the element of free will present in all human 
decisions by which we do not follow blindly the opinion dictated by our relationships. 
Our aim is to analyze how the interplay between this free will and the interactions 
amongst individuals affects group formation in opinion dynamics. In the language of 
statistical mechanics, what we are doing is to add noise to the deterministic dynamics 
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and analyze which aspects of the model are robust against the introduction of noise. 
Noise is introduced by allowing an individual opinion to change to another randomly 
chosen value in the whole opinion space. Under some circumstances this turns out to be 
equivalent to allowing each agent to return, at some random times, to a specific opinion 
preferred by him. 

Our analysis, based upon numerical integrations of the corresponding master equa- 
tion as well as Monte-Carlo simulations, reveals new and interesting phenomenology. 
There exists a critical value rric of the noise intensity, which depends on the confidence 
range, such that for noise larger than this value the system becomes disorganized and 
group formation does not occur. We provide a linear stability analysis that reproduces 
the order-disorder transition that occurs at rric. For noise smaller than mc, the steady- 
state probability distributions in opinion space broaden with respect to the noiseless 
case, but still have large peaks and group formation can be unambiguously defined by 
looking at the maxima of the distributions. The group formation occurs by a series of 
bifurcations that mimic those that occur in the noiseless case. 

An important aspect of our work, that we want to stress here, is that the numerical 
Monte-Carlo simulations do not necessarily agree with the results of the master equation. 
This is due to the inherent finite-size-induced fiuctuations that occur in the simulations. 
A similar warning is required when one tries to infer about possible applications of the 
model to real situations. For example, it is possible to find regions of bistability where 
dynamical transitions between a single group and two groups occur. These transitions 
do not occur in the infinite-size thermodynamic limit taken routinely in most studies. 
This stresses the role that a finite size has on the dynamics of social systems f25]. 

This paper is organized as follows. Deffuant et al. model is briefiy reviewed in 
Sec. [2} The main results are presented in section [3| devoted to study this model in the 
presence of noise. In Sec. [4] we use a linear stability analysis to derive approximately 
the critical value of the noise intensity for the formation of opinion groups. Summary 
and conclusions are presented in Sec. [5| 

2. Review of Deffuant et al. model 

Let us consider a population with individuals. We will denote by the number 
representing the opinion on a given topic that individual i has at time-step n. As 
mentioned in the introduction, the opinion is a real variable in a finite interval and, 
without loss of generality, we take G [0, 1]. Initially, it is assumed that the values Xq 
for i = 1, . . . , A are randomly distributed in the interval [0, 1]. Dynamics is introduced 
to refiect that individuals interact, discuss, and modify their opinions. In the original 
version of the model [6j, at time-step n two individuals, say i and j, are randomly 
chosen. If their opinions satisfy — < e, so that they are close enough, they are 
modified as: 



x^ 

j 
x^ 



X^ + 




(1) 
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otherwise they remain unchanged. Whether the opinions have been updated or not, 
time increases n ^ n + 1. As a consequence of the iteration of this dynamical rule, 
the system reaches a static final configuration which, depending on the values of the 
parameters e and /i, can be a state of full consensus where all individuals share the same 
opinion, or of fragmentation with several opinion groups. It is customary to introduce 
the time variable t = nAt, where At = measuring the number of opinion updates 
per individual, or number of Monte-Carlo steps (MCS). 

The parameter fi is restricted to the interval (0, 0.5]. It determines the convergence 
time between individuals as well as the number of final groups [23l [24] . For small values 
of /X, the individuals slightly change their opinions during meeting, while for /x = 0.5 the 
interacting individuals fully compromise and, after the meeting, they share the same 
opinion. As in most studies, we will adopt from now on in this paper the value /x = 0.5. 
The parameter e, which runs from to 1, is the confidence parameter. Starting from 
uniformly distributed random values for the initial opinions, the typical realization is 
that for large values, e > 0.5, the system evolves to a state of consensus where all 
individuals share the same opinion and that, decreasing e, the population splits into 
opinion groups separated by distances larger than e. 

The process can be described in terms of a master equation for the probability 
density function for an individual opinion x at time t. The Appendix contains a 

derivation of this master equation in the presence of the additional noise term described 



in section |3j Equation (A.4) with m = is the master equation for the noiseless 
original Deffuant et al. model, first obtained in fT2]. A detailed analysis based on 
its numerical integration [12[ |20|f] shows that there are four basic modes of group 
appearance, dominance, or splitting, which are called bifurcations ^[[20] in this context 
(see Fig.jl]): nucleation of two minor groups symmetrically from the center of the opinion 
interval (type-1, such as the birth of two minor groups from the boundaries at e = 0.5 
in Fig. [T]); nucleation of two major groups from the central one (type-2, as occurring 
at e ^ 0.266 in Fig. [T]); nucleation of a minor central group (type-3, as occurring at 
e ^ 0.222 in the figure); and, finally, sudden increase of the mass of that central group 
accompanied by a sudden drift outwards in the location of the two major groups (type- 
4, at e 0.182 in the figure). In this sequence "major" opinion groups contain a high 
fraction of the population, while "minor" groups contain a much smaller fraction (of 
the order of 10~^ or smaller). The bifurcation pattern repeats itself as e decreases even 
further. 

It is important to emphasize that the situation depicted in Fig. [T] is the result 
for steady solutions of the master equation attained at long times starting from a 
uniform initial distribution. Many other steady solutions of the master equation exist. 
In particular, any combination of delta-functions is a steady solution of the noiseless 
master equation provided they are separated by more than a distance e. We stress 
also that this analysis based upon the master equation corresponds to the limit case 

if In order to compare with the results of reference [12 we note that our parameter e is related to their 
parameter A by e = 1/2A. 
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where the number of individuals tends to infinity. In Monte-Carlo simulations of 
the microscopic rules, or in practical applications with a necessarily finite number of 
individuals, some features need to be considered. It is still true that each realization 
ends up in a small number of groups (both of major and minor type), all individuals 
within a group holding exactly the same opinion. However, the exact location of the 
groups might vary with respect to the master equation prediction and minor groups 
might not appear depending on the particular realization and the total size of the 
population. Furthermore, there could be realizations in which even the number of 
observed major groups differs from the prediction of the master equation. These effects 
are more pronounced the smaller the number of individuals. In the same Fig. [T] we 
have plotted the distribution of observed groups, averaged over many realizations for 
two different number of individuals A^, where the aforementioned properties can clearly 
be observed. For instance, for e = 0.28, the master equation predicts that there should 
be only one major group, centered at x = 0.5, and two minor groups. However, in 
almost half of the realizations with = 1000 individuals the opinions split instead in 
two major groups centered around x = 0.28 and x = 0.72 and, eventually, some minor 
extreme groups. 



3. Effect of noise 



Noise is introduced as a random change of an individual's opinion. Specifically, we 
modify the dynamics as follows: at time-step n the original dynamical rule, Eq. ([T]), 
applies only with probability 1 — m. Otherwise, a randomly chosen individual i changes 
his opinion to a new value xl^^-^ drawn from a uniform distribution in the interval [0, 1] 
and all other opinions remain unchanged. The probability m is a measure of the noise- 
intensity. Note that, quite generally, this rule is equivalent to allowing each agent to 
return to a specific, basal, opinion preferred by him, provided that the basal opinions 
are randomly distributed amongst the agents. In this section, we will study in detail 
the effect of this new ingredient in the dynamical evolution of the model. We analyze 
both the results coming from a numerical integration of the master equation as well as 
numerical Monte-Carlo-type simulations of the microscopic rules of the model. 



3.1. Master equation approach 



The Appendix contains a derivation of the master equation, Eq. (A.4), appropriate for 
this process. We have first obtained the asymptotic distribution Poo(^) = lim^^oo P{x^ t) 
of the master equation starting from a suitable initial condition t = 0). As in other 
studies, we assume that the initial condition represents a uniform distribution in opinion 
space, i.e. P{x^t = 0) = 1 for x G [0, 1] and P{x^t = 0) = otherwise. For m = 
the steady-state distribution Poo(^) it is a sum of delta- functions located at particular 
points. In the case m > the steady distributions are no longer delta-functions but 
still are peaked around some particular values if e is not too small or m not too large. 
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Figure 1. Top panel: bifurcation diagram of the noiseless Deffuant et al.^ model 
reproduced from reference [22], coming from a numerical analysis of the master 
equation of the model J2l [21]. Lines show the position of the opinion groups as a 
function of e. By arrows we indicate the location of the four basic types of bifurcations: 
type-1 at e ~ 0.5, type-2 at e ~ 0.266, type-3 at e ~ 0.22, type-4 at e ~ 0.182, and this 
pattern repeating as e decreases even further [21 . In the bottom panels we plot in a 
logarithmic grey-scale the asymptotic probability distribution Poo{x) (values smaller 
than 2 x 10~^ are plotted white), as a function of e, resulting from extensive numerical 
simulations of the microscopic model with N = 10^ and N = 10^ agents, respectively. 
Note that a single realization using the Monte-Carlo microscopic rules of the model 
leads to a probability distribution which is a sum of delta-functions. The distributions 
displayed in the panels are the result of an average over 10^ realizations for = 10^ 
and 2 x 10^ realizations for N = 10^ and a histogram bin size Ax = 0.01. In all cases, 
both in the master equation as in the Monte-Carlo simulations, the initial condition 
represents opinions which are random and uniformly distributed in the interval [0, 1]. 
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Figure 2. Plot, in a logarithmic grey scale, of the asymptotic probability distributions 
Poo(^) as functions of e, obtained after a numerical integration of the master equation 
( |A.4[ ) for m = 0.01, top panel, and m = 0.1, bottom panel. We have used a fourth-order 
Runge-Kutta method with a time step At = 0.1 for the time evolution and Simpson's 
rule for the integrals in x-space with a discretization Ax = 1/M, M = 2000. We have 
checked in some tests cases that smaller time or discretization steps do not change 
significantly the results. As in Fig {l] we show by arrows the location of the bifurcation 
points. 



We have plotted in Fig. [2] the master equation steady probabihty distributions Poo{x) 
as a function of the parameter e for two different values of the noise intensity. In 
the small noise case m = 0.01 it is still possible, for not too small e, to identify the 
same type of bifurcations than in the noiseless case by looking at the maxima of the 
probability distributions: a type-1 bifurcation at e = 0.5 where minor groups begin to 
form at X = and x = 1; a type-2 bifurcation at e 0.2695 where the distribution 
switches from having one single maximum at x = 0.5 to having two maxima of equal 
height located at x 1/4 and x 3/4; a type-3 bifurcation at e 0.250 where a 
central maximum begins to grow; a type-4 bifurcation at e 0.1835 where three equally 
spaced maxima of equal height atx?=^l/6, x = l/2 and x 5/6 appear. This pattern of 
bifurcations repeats as e decreases even further. However, type-1 and type-3 bifurcations 
are somewhat ambiguous to define since the relative importance of the minor groups 
actually increases continuously instead of sharply increasing when new maxima begin to 
form. It is worth stressing that, at variance with the noiseless case, the location of the 
major groups (defined as the absolute maxima of the distribution) does not vary with 
e until a new bifurcation of type-2 or type-4 is reached. These maxima are regularly 
located at x l//c, 3//c, . . . , (A; — l)//c for = 2, 4, 6, . . . 

The same general structure can be observed in the case of larger noise m = 0.1, 
although the distributions are much broader now. The location of the main bifurcation 
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points can be located at e 0.2965 (type-2) and e ^ 0.1755 (type-4). The type-1 and 
type-3 transitions are very imprecisely defined, specially for smaller values of e. 

For both noise values, one observes that groups become less defined and finally 
are replaced by a more or less unstructured distribution for e below a critical value 
Cc which increases with m. Alternatively, one realizes the existence of a critical value 
rric = mc{e)^ increasing with e, above which the group structure disappears from the 
steady distribution. 

A somehow expected feature that emerges from the data shown in Fig. [2] is that the 
width of the steady distributions grows with the noise intensity. An explicit expression 
for the width of the single maximum present when e > 1 (so that all individuals are 
allowed to interact) could be obtained from the master equation, since in this case the 
moments form a closed hierarchy. Defining the moments Mi and M2 as in the Appendix, 
the variance = M2 — satisfies: 

In this limiting case (where the main feature of the model, bounded confidence, has 
been lost since everybody is able to interact with everybody), the variance reaches a 
steady-state in which the width increases with noise as a ~ m^/^. 



3.2. Comparison with Monte- Carlo simulations 

Once the master equation predictions have been established, and before comparing with 
the results coming from the Monte-Carlo numerical simulations using the microscopic 
rules, a word of warning is, as in the noiseless case, required. In g simulations with a 
finite number of individuals the dynamics of the probability distribution as well as 
its asymptotic, steady-state, values might not coincide with the analysis of the master 
equation. We have found this deviation to be more pronounced in the case of being 
close to a bifurcation point. For example, in Fig. [3] we plot the time evolution of the 
probability coming from Monte-Carlo simulations of the model for diflFerent system sizes 
and the results of the master equation in the case e = 0.28, close to a type-2 bifurcation 
point. It can be seen that, although the Monte-Carlo simulation and the master equation 
agree initially very well, they start to deviate after a time that depends on the number 
of individuals N: the larger A, the longer the time than the Monte-Carlo simulations 
are faithfully described by the master equation. In this particular case, e = 0.28, it can 
be seen that the master equation agrees with the Monte-Carlo simulations up to a time 
t ^ 10 for A = 10^ and a time t ^ 100 for A = 10^. In view of this difference, it is 
surprising that the Monte Carlo steady-state distributions show only small (although 
observable) finite-size effects. Quite similar functions describe the steady-state data for 
both A = 10^ and A = 10^, see Fig. HI As can be seen in Fig. [sj however, while 
the numerical solution of the master equation tends to the steady-state distribution 
^00(^)5 the Monte-Carlo simulations tend to another distribution, Pst{x). These two 
distributions are very different: Poo{x) has a large maximum (large group) at x = 0.5 
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Figure 3. Probability distribution function P{x^t)^ for intermediate time steps, from 
Monte-Carlo (MC) simulations (histograms binned with bin size Ax = 5 x 10""^) for 
two different system sizes N = 10^ and N = 10^, and the master equation (ME) 
integrations of Eq. (A.4) starting with a flat distribution with e = 0.28 and m = 0.01. 
The distributions are an average over 10^ realizations for N = 10^ and 10^ realizations 
for TV = 10^ 



and two much smaller maxima x ^ 0.127 and x ^ 0.873, whereas Pst{x) has two 
equal maxima at x 0.25 and x ^ 0.75. Although surprising at first, it turns out that 
the steady-state distribution Pst{x) coming from the Monte-Carlo simulations of the 
model is also very close to a steady-state solution of the master equation (A.4) having 
two major groups. However, this last steady-state solution can not be obtained as an 
asymptotic solution of the master equation lim^^oo P{x^t) starting from a uniform initial 
condition P{x^t = 0) = 1 for x G [0, 1]. It turns out that it is reached when starting 
instead from an initial condition asymmetric with respect to the center of the interval. 

Summing up, for e = 0.28 there are two steady-state solutions of the master 
equation, Pst{x) and Poo(^)- Starting from a uniform initial condition, Poo(^) is the 
one asymptotically reached as a solution of the master equation. However, Pst{x) is, 
up to finite-size effects, the one reached instead in the Monte-Carlo simulations. We 
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Figure 4. The symbols denote the steady-state probabihty distribution function 
Pstix)^ obtained from Monte-Carlo simulations for system sizes N = 10^ (•) and 
N = 10^ (x) in the case e = 0.28, m = 0.01. The dashed line is obtained as 
an asymptotic solution of the master equation starting with a non-uniform initial 
condition. These three distributions are very similar in this logarithmic scale, although 
there are differences in the height of the maxima (see Fig. [3|. It turns out that 
close to each of the two maxima, the distributions can be well fitted by a Lorentzian 
function [19 . The solid line is the distribution Poo{x) coming from a numerical solution 
of the master equation Eq. (A.4) starting with a uniform distribution. 



interpret this in terms of the relative stability of both solutions: introducing a fluctuation 
5P{x) on the solution Poo(^) it is then possible to reach the solution Pst{x)^ but not 
the reverse. This fluctuation 5P{x) needs to be asymmetric, 5P{x) ^ bP{\ — x) and it 
appears naturally in Monte-Carlo simulations because of the finite number of individuals 
and it is more probable the smaller the value of A^. This explains that the system 
with smaller deviates earlier from the solution of the master equation. If one induces 
artificiallj|§] such a non-symmetric perturbation in the solution Poo(^) or, alternatively, 
one starts with a non-uniform, asymmetric initial condition, then the master equation 
tends to Pst[x). 

The existence of more that one steady solution of the master equation seems to 
be a general feature. In fact, if the distributions in Fig. [2] are recalculated by slowly 
increasing and decreasing e without resetting the initial condition to = 0) = 1 

after each change in e, we observe the hysteresis behavior typical of bistability occurrence 
close to first-order transitions. Which one of the possible steady solutions is observed in 

§ If one is not careful enough, the perturbation might also appear as a numerical instability of the 
integration method. 



Noisy continuous-opinion dynamics 



11 



10.0 



1.0 



0.1 



I 



1.0 
0.8 

0.6 t 



0.2 


















e = 0.28 



0.2 
0.0 

1.0 
0.8 

0.6 

0.4 

0.2 
0.0 

1.0 
0.8 

0.6 

0.4 

0.2 
0.0 



G = 0.31 



6 = 0.35 



e = 0.80 







200000 400000 600000 

t (MCS) 



800000 



Figure 5. Time series of the opinion distributions P(x, t) at four values of e and 
m = 0.01 and N = 1000. Grey scale is logarithmic, with values smaller than 0.1 
plotted in white. Simulations are initialized with random homogeneous distributions 
1000 MCS before the first time shown in the plots. The three first panels are near a 
type-2 bifurcation. At e = 0.28 the system appears polarized in two opinion groups. 
At e = 0.35 a single major opinion has been established, with two lateral minor groups. 
At e = 0.31 the systems fluctuates between these two states. A single group is observed 
at large e (bottom panel). 



the Monte-Carlo simulations depends on the parameter e. It could even happen that the 
inherent fluctuations of a finite system take it from one solution to another and back. 
This sort of bistability is observed, for instance, in Monte-Carlo simulations at e = 0.31. 
As shown in Fig. [5| the evolution displays multiple jumps between two solutions: one 
with two maxima of equal height and another one with a large central maximum and 
two smaller maxima near the edges of the opinion interval. 

Looking at Fig.[2| one can see that for small e or, alternatively, for noise intensity m 
larger than a critical value rric which increases with e, the bifurcations become blurred 
and the maxima of the distributions are not evident, implying the inhibition of group 
formation. This happens at e ^ 0.1 for m = 0.1 and at e ^ 0.01 for m = 0.01. A 
similar effect can be observed in the Monte-Carlo simulations and can be described in 
terms of an order-disorder transition: order identified with the state with well defined 
opinion groups and disorder identified with the state without groups. To define in 
a more quantitative way this transition, we have used the so-called group coefficient 
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which aims to characterize the existence of groups [26]. The definition of the 
group coefficient Gm starts by dividing the opinion space [0, 1] in M equal boxes and 
counting the number of individuals li which, at time step n, have their opinion in the box 
[{i — 1)/M, i/M]. One next introduces an entropy-like measure Sm = ~ X^^^i n^^n- 
Note that < Sm < InM, and that the minimum value Sm is obtained when all the 
individuals are in just one box, while the maximum value, Sm = InM, is reached when 
li = N/M^ i.e. when the opinions are uniformly distributed in the interval [0, 1]. Finally, 
the opinion group coefficient is defined as [26] 

GM = M-'(^e^^y (3) 

where the over-bar denotes a temporal average in steady conditions and (•) indicates 
an average over different realizations of the dynamics. Note that 1/M < Gm ^ 1- 
Large values Gm ^ 1 indicate that the opinions are evenly distributed along the full 
opinion space (a situation identified with disorder), while small values of Gm indicate 
that opinions peak around a finite set of major opinion groups (a situation identified 
with order). 




Figure 6. Opinion group coefficient Gm=ioo versus m for e = 0.05, 0.2, 0.4, and 
0.6, from left to right, respectively as obtained from Monte-Carlo simulations with 
N = 10^ individuals (dots, the solid line is a guide to the eye). The location of the 
order-disorder transition is defined as the value rric for which Gm =0.9. Note that rUc 
is an increasing function of e and that the transition is very abrupt for small e. 

The data coming from Monte-Carlo simulations, see Fig. [6j show that Gm is an 
increasing function of the noise intensity m and saturates to its maximum value Gm ^ 1 
for large enough values of m. The transition from group formation to disorder will 
be defined, somehow arbitrarily but precisely, as the value rric of the noise intensity 
for which the group coefficient reaches the value Gm = 0.9. For small values of the 
confidence parameter e the transition to the homogeneous state is abrupt and occurs for 
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small values of m. If one increases e, the transition becomes less abrupt and a higher 
noise intensity is needed to obtain the homogeneous, group-free, state. This last feature 
can be explained by the linear stability analysis that we shall develop in the next section. 



4. Linear stability analysis 

We have shown that opinion groups still form in presence of small amounts of noise, but 
an unstructured state without groups dominates the opinion space for noise larger than 
a critical noise intensity rric. Although the transition to group formation is a nonlinear 
process, one can still derive approximate analytical conditions for the existence of group 
formation in the parameter space (e, m) by performing a linear stability analysis of the 
unstructured solution of Eq. ( A. 41 ). This is greatly simplified if one neglects the infiuence 



of the boundaries and assumes that the interval [0, 1] is wrapped on a circle, i.e. there are 
periodic boundary conditions at the ends of the interval. This would be a reasonable 
approximation to describe the distribution far from the boundaries if e, which fixes 
the interaction range, is sufficiently small. In this case the homogeneous configuration 
Ph{x) = 1 is an approximation to the unstructured steady solution of the master 
equation. Analysis of its stability begins by introducing = Ph{x) + ^4^6^^^+^^^, 

where q is the wave number of the perturbation, its growth rate and Aq the amplitude. 



After introducing this ansatz in Eq. (A.4) we find the dispersion relation giving the 
growth rate of mode q: 



Xq = 4e(l - m) 



4sin(ge/2) sin(ge) 



m. (4) 



qe qe 

This is plotted in Fig. [7](a), for several values of m. The maximum value of occurs at 
Qmax = 2.7907/6. It turns out that he maximum value Xq^^^ is negative for m > rric and 

positive for m < rric = with a ^ 0.8676. Alternatively, for fixed m the maximum 

a + e 

avfi 

growth rate is negative for e < ec = , and positive for e > ec. Therefore, the 

1 — m 

homogeneous state is unstable and group formation is possible only for m < rric or 
e > ec- The numerical values are ec = 0.096 for m = 0.1 and 6c = 0.0088 for m = 0.01, 
in reasonable agreement with the behavior observed for the master equation dynamics 
in Fig. [2| Comparison with Monte-Carlo simulations is performed in Fig. [7](b) where we 
plot the critical value rUc obtained from the group coefficient Gm as described earlier. 
We see in the figure that the agreement is very good for small e but deviates for larger 
values. This is consistent with the fact that neglecting boundary effects is expected to 
be appropriate only for small e. Finally, it is possible to estimate roughly the number of 
groups n by a simple argument: n is related to the wavelength of the maximum growth 
as n = qmax/'^TT OT n = 0.444/e for e > ec- This result is in qualitative agreement with 
the ^-rule, which says that the number of major groups after group formation is roughly 
determined as the integer part of ^ (see [20l[21] for details). 
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Figure 7. (a) Growth rate, Eq. Q, of a perturbation to the homogeneous solution as 
a function of qe for e = 0.35 and m = 0.0, 0.10, 0.287, and 0.35 from top to bottom. It 
shows that Xq becomes negative for high values of m. (b) Phase diagram on the plane 
(e, m) as coming from Monte-Carlo simulations with N = 10^ individuals (dots), using 
the criterion based on the group coefficient Gm=ioo as described in the text. The solid 
line is the prediction of the linear stability analysis. Groups appear below the lines, 
whereas the unstructured state is stable above. 



5. Summary and conclusions 

In this paper we have studied the Deffuant et al. model for continuous-opinion dynamics 
under the presence of noise. Besides the usual rules of the model, we give each individual 
the opportunity to change, with a certain probability m, his opinion to a randomly 
selected opinion inside the whole opinion space. The final behavior depends of the 
confidence or interaction parameter e and the noise intensity m. 

We have first reviewed the original noiseless version of the model. We have shown 
that for small number of individuals and depending of the particular realizations, 
the exact location of the opinion groups might vary with respect to the predictions of 
the master equation. In particular, minor groups might not appear and there could 
be realizations in which even the number of observed major groups differs from the 
prediction of the master equation. 

We have derived (Appendix) a master equation for the probability density function 
t) which determines the individuals density or distribution in the opinion space. 
Numerical integration of this equation from uniform initial conditions reveals that for 
m > the steady distributions are no longer delta- functions as for m = 0, but are still 
peaked around some well defined maximum value, with some non-vanishing width. By 
looking at those maxima we are able to identify the same type of bifurcations than in the 
noiseless case |lT2l [20] . At variance with the noiseless case, the location of the maxima 
(the central opinion of the groups) does not depend on e until a new bifurcation point 
is reached. 

We have also found that in the noisy case the asymptotic steady-state probability 
distributions reached by Monte-Carlo simulations might not coincide with the ones 
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obtained from the master equation dynamics starting from the same symmetric initial 
condition. This deviation is more pronounced in the case of being close to a bifurcation 
point. In particular, we have presented a situation where, starting from a uniform 
initial condition, a particular stationary distribution, Pqq{x)^ is actually reached by the 
master equation but another distribution , Pst{x)^ is the one reached instead in Monte- 
Carlo simulations. The time in which the Monte-Carlo simulations begin to deviate 
from the master equation depends on system size: the smaller the size, the earlier the 
deviation occurs although the final Monte-Carlo distribution Pst{x) shows only small 
size effects. Remarkably, Pst{x) turns out to be close to another steady solution of the 
master equation. Thus, the discrepancy observed during the dynamics does not seem 
to be simply a trivial finite- effect. We interpret it in terms of the relative stability of 
both solutions by adding an asymmetric perturbation to Poo{x). It is then possible to 
reach the solution Pst{x) obtained from Monte-Carlo simulations, but not the other way 
around . Asymmetric fluctuations appear naturally in Monte-Carlo dynamics because 
of the finite number of individuals, and are larger for smaller A^. We have also shown 
that the fiuctuations present in Monte-Carlo simulations are even able to induce jumps 
from one solution to another and back. 

An order-disorder transition to group formation induced by noise has been 
characterized using the so-called group coefficient Gm for the simulations performed 
with Monte-Carlo dynamics. We have found that Gm is an increasing function of the 
noise intensity m and saturates to its maximum value Gm ^ 1 for large enough values of 
m. For small values of the confidence parameter e the transition to the disordered state 
is abrupt and occurs for small values of m. If one increases e, the transition becomes 
less abrupt and a higher noise intensity is needed to obtain this state. Using a linear 
stability analysis of the unstructured (no groups) solution of the master equation we 
have derived approximate conditions for opinion group formation as a function of the 
relevant parameters of the system. We have found qualitative agreement between the 
linear stability analysis and numerical simulations. The agreement is better for small 
values of e where boundary effects, neglected to make feasible the linear analysis, are 
less important. However, we should emphasize that the pattern selection of this model 
is, with noise and without it, intrinsically a nonlinear phenomenon and obtaining the 
exact critical conditions for opinion group formation remains a challenge. 

Our work stresses the importance that fiuctuations and finite-size effects have 
in the dynamics of social systems for which the thermodynamic limit is not justified 
[25]. Further work will address the effect that these ingredients have in the dynamics 
of continuous-opinion models in the presence of an external infiuence, or forcing, 
representing the role of advertising. 
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Appendix 



Here we derive the master equation, i.e. the evolution equation for Pn{x)^ the probabihty 
density function (pdf) of the opinions at step n for the model introduced in this paper. 
Note that Pn{x) is constructed from the histogram of all individual opinions at step 
n. Let us first find the evolution of the pdf for those two particular individuals i^j 
that have been selected for updating at step n according to the basic rule Eq. ([T]). We 
will denote by Pn+i{x) the pdf of the opinion of individual i at the step n + 1, i.e the 
probability that x^^^ adopts the value x. According to that rule, it is straightforward 
to derive the evolution equation 



dx ^dx'^Pji {x^Pji {x'!^5 



X 



rpl _|_ np'J 



+ / dx^dx^^Pji(^x^)Pji(^x^^)S(^x x^)^ (^-l) 

and a similar expression for P^^-^{x). The integrals over and x^n run both over the 
interval [0, 1]. In this equation an independence approximation for the variables 
has been implicitly assumed, i.e. their joint pdf is supposed to factorize as 

Pn{x\^xi^) = Pn{x\)Pn{xi^). TMs is au uncontrolled approximation whose validity can 
only be established by an ulterior comparison with the Monte-Carlo simulation of the 
microscopic rules. For the individuals k ^ i^j whose opinion does not change at time 
step n + 1 we have simply P^^-^{x) = P^{x)- 

The pdf Pn+i{x) has several contributions: (i) With probability 1 — m two 
individuals, say i, j, are chosen for updating according to the basic evolution rule Eq. ([T]) 
and N — 2 variables remain unchanged, (ii) With probability m one individual is chosen 
for updating according to the noise rule and N — 1 variables remain unchanged; the new 
opinion of the selected individual is sampled from an, in principle, arbitrary distribution 
Pa{x)^ although in this paper we have taken throughout that Pa{x) is the uniform 
distribution Ph{x) = 1 in the interval [0, 1]. After consideration of these contributions 
we are led to the evolution equation 



(x) 



m] 



N-2 



N 



+ m 



N-1 
N 



Pn{x) + -P, 



[X 



(A.2) 



Replacing and P^^i{x) from Eq. (A.l) one obtains after some algebra 

(1 — m) 



Pn+lix) = Pnix) + 
-2Pn{x) 



N 



4 / dx'Pn{2x - x')Pn{x') 

\x-x'\<e/2 



dx'Pnix') 



\x—x'\<e 



m 



+ - [Pa{x) - P^{X)] . 



(A.3) 
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: / dx'P{2x-x',t)P{x',t) 

J\x-x'\<e/2 

m[Pa{x)-P{x,t)], (A.4) 



The integrals over x^ run over the interval [0, 1], and it has to be imposed that Pn{x) = 
if X ^ [0, 1]. We now take the continuum limit Pn{x) P{x^t) with a time t = nAt 
and taking the limit At = ^ as ^ oc, to obtain: 

dP{x,t) 

= — III ) 4: 

^\x-x'\<e/2 

-2P{x,t) I dx'P{x',t) 

J\x-x'\<e 

which is the master equation of the Deffuant et al. model in the presence of noise and 
the basis of our analysis. The noiseless case, m = 0, was first obtained in reference \TT\. 
We note here the symmetry property of the master equation: if the initial condition is 
symmetric around the central point x = 1/2, namely that P(x,t = 0) = P(1 — x,t = 0), 
then this property holds for any later time, = P{1 — x^t)yt > 0. 

The time evolution of the first moments of t) can be computed from the master 

equation. Defining the moments as Mk{t) = J dxx^P{x, t) one finds easily that = 

(normalization condition) and that the first moment evolves as = ^ {^i ~ ^i)^ 

being the first moment of the distribution Pa{x). Therefore, if m > 0, the average 
opinion tends to Mi = independently of the initial condition, and it is always 
conserved in the noiseless case m = 0. Expressions for higher-order moments can only 
be obtained in the special case e > 1, as discussed in the main text. 
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